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Abstract: 

The dominant production process for heavy charged Higgs bosons at the LHC is the associated 
production with heavy quarks. We have calculated the next-to-leading-order supersymmetric 
QCD corrections to charged-Higgs production through the parton processes qq, gg — > tbH^ 
and present results for total cross sections and differential distributions. The QCD corrections 
reduce the renormalization and factorization scale dependence and thus stabilize the theoretical 
predictions. We present a comparison of the next-to-leading-order results for the inclusive cross 
section with a calculation based on bottom-gluon fusion gb tH^ and discuss the impact of 
the next-to- leading-order corrections on charged-Higgs searches at the LHC. 



1 Introduction 



The Higgs mechanism [1] is a cornerstone of the Standard Model (SM) and its supersymmetric 
extensions. The masses of the fundamental particles, electroweak gauge bosons, leptons, and 
quarks, are generated by interactions with Higgs fields. The search for Higgs bosons is thus 
one of the most important tasks for high-energy physics and is being pursued at the upgraded 
proton-antiproton collider Tevatron with a centre-of-mass (CM) energy of 1.96 TeV, followed by 
the proton-proton collider LHC with 14 TeV CM energy scheduled to start taking data in 2010. 

The minimal supersymmetric extension of the Standard Model (MSSM) requires two Higgs 
doublets leading to five physical scalar Higgs bosons: two (mass-degenerate) charged particles 
H±, one CP-odd neutral particle A, and two CP-even neutral particles h,H. The discovery of a 
charged Higgs boson, in particular, would provide unambiguous evidence for an extended Higgs 
sector beyond the Standard Model. Searches at LEP have set a limit Mjj± > 79.3 GeV on the 
mass of a charged Higgs boson in a general two-Higgs-doublet model [2]. Within the MSSM, 
the charged-Higgs mass is constrained by the pseudoscalar Higgs mass and the W-boson mass 
through M^± = -|- at tree level, with only moderate higher-order corrections [3-6]. 
A mass limit on the MSSM charged Higgs boson can thus be derived from the limit on the 
pseudoscalar Higgs boson. Ma > 93.4 GeV [7], resulting in Mh± > 120 GeV. At the Tevatron, 
searches for light charged Higgs bosons in top-quark decays t — > bH='= [8, 9] have placed some 
constraints on the MSSM parameter space, but do not provide any further generic bounds on 
M±. 

The LHC will extend the search for charged Higgs bosons to masses up to Afjji ^ 600 GcV [10, 
11], where the reach depends in detail on the values of the supersymmetric parameters. In this 
paper we shall focus on the most promising search channel for heavy H='= (with Mh± ^ mt) at 
the LHC, which is the associated production of charged Higgs with heavy quarks, 

pp^tbH±+X. (1.1) 

Alternative production mechanisms like quark-antiquark annihilation qq' — > H^, H^ + jet pro- 
duction, associated H='=W^ production or Higgs pair production have suppressed rates, and it 
is not yet clear whether a signal could be established in any of those channels (see Ref. [12] 
and references therein). Some of the above production processes may, however, be enhanced in 
models with non-minimal flavour violation (see, e.g., Ref. [13]). 

Two different formalisms can be employed to calculate the cross section for associated tbH='= 
production. In a four-flavour scheme (4FS) with no b quarks in the initial state, the lowest-order 
QCD production processes are gluon-gluon fusion and quark-antiquark annihilation, gg tbH^ 
and qq tbH^, respectively. The inclusive cross section for gg — > tbH^ develops potentially 
large logarithms oc ln(/[/F/mb), which arise from the splitting of incoming gluons into nearly 
collinear bb pairs. The large scale /Up of (9(Mh±) corresponds to the upper limit of the collinear 
region up to which factorization is valid. The ln(;UF/'"T^b) terms can be summed to all orders 
in perturbation theory by introducing bottom parton densities. This defines the so-called five- 
flavour scheme (5FS) [14]. The use of bottom distribution functions is based on the approximation 
that the outgoing b quark is at small transverse momentum and massless, and the virtual b 
quark is quasi on-shell. In this scheme, the leading-order (LO) process for the inclusive tbH^ 
cross section is gluon-bottom fusion, gb tH^. The next-to-leading order (NLO) cross section 
in the 5FS includes 0{as) corrections to gb — s- tH^ and the tree-level processes gg — > tbH^ and 
tbH=^. 

To all orders in perturbation theory the four- and five-flavour schemes arc identical, but the 
way of ordering the perturbative expansion is different, and the results do not match exactly at 
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finite order. For the inclusive production of neutral Higgs bosons with bottom quarks, pp — > 
bbH + X, the four- and five- flavour scheme calculations numerically agree within their respective 
uncertainties, once higher-oder QCD corrections are taken into account [15-18]. However, no 
NLO comparison of the 4FS and 5FS calculations for charged-Higgs production with heavy 
quarks exists so far. 

There has been considerable progress recently in improving the cross-section predictions for 
the associated production of charged Higgs bosons with heavy quarks by calculating NLO-QCD 
and SUSY-QCD corrections in the four and five-flavour schemes [19-24]. The inclusion of higher- 
order effects is crucial for an accurate theoretical prediction and, eventually, a determination 
of Higgs-boson parameters from the comparison of theory and experiment. In this paper we 
present an independent calculation of the NLO supersymmetric QCD corrections to the process 
pp — > tbH^ + X in the 4FS. The calculation within the 4FS allows to describe the dynamics of 
the final-state bottom quark, which in the 5FS scheme calculation at LO is assumed to be always 
produced at small transverse momentum and is thus treated inclusivelj0. However, Monte Carlo 
simulations show that in about 20% of pp — > tbH^ + X events at the LHC the b quark from 
the production process has a transverse momentum larger than the b quark from the top-quark 
decay, and will thus contaminate the event reconstruction [25]. We therefore provide state-of- 
the art NLO predictions not only for the inclusive cross section but also for various differential 
distributions. In contrast to previous analyses our results are based on the consistent use of a 
four-flavour parton distribution function. Furthermore, we present the first comparison of the 
4FS and 5FS calculations at NLO for the inclusive tH^ cross section. 

The paper is organized as follows: In Section [2] we shall describe the calculation of the NLO 
supersymmetric QCD corrections. Numerical results for MSSM Higgs-boson production at the 
LHC are presented in Section [3l We conclude in Section [H The Appendix provides details on 
the scenario of the supersymmetric model under consideration. 

2 Calculation 

2.1 LO processes and conventions 

In the 4FS the production of charged Higgs bosons in association with top and bottom quarks 
proceeds at LO through the parton processes [26-28] 

gg tbH~ and qq tbH" , (2.1) 

and the charge-conjugate processes with the tbH"*" final state. Throughout this paper we present 
results for the tbH~ channels, unless stated otherwise. Generic Feynman diagrams that con- 
tribute to the LO processes (|2.ip are displayed in Fig. HJ^a). 

In the MSSM, the Yukawa coupling of the charged Higgs boson H~ to a top and bottom 
quark is given by 

5tbH- = ^ {"^Pr cot (3 + tan /?) , (2.2) 

where v = ■\/v'l + V2 = (^/2Gf)~^/2 is the vacuum expectation value of the Higgs field in the 
Standard Model, and Gp = 1.16637 x 10^^ GeV^^ [29] is the Fermi constant. The ratio of the 
vacuum expectation values vi and V2 of the two Higgs doublets is denoted by tan (5 = V2/V1, and 
^L/i?, = (1 T 75)72 are the chirality projectors. 



^This shortcoming of the 5FS, however, is rectified when going to NLO, where the process gg — > tbH^ con- 
tributes as part of the real corrections. 



2 




Figure 1: A generic set of diagrams (a) for the Born level, (b) for virtual gluon exchange, (c) 
virtual gluino and squark exchange, (d) gluon radiation, and (e) gluon-(anti) quark scattering in 
the subprocesses qq, gg tbH~, etc. 



3 



2.2 NLO supersymmetric QCD corrections 

The NLO supersymmetric QCD corrections comprise virtual one-loop diagrams, Fig. mb,c), 
gluon radiation processes, Fig. [^d), and gluon-(anti)quark scattering reactions. Fig. [T]^e). The 
NLO QCD calculation of the SM processes qq, gg — > QQH, where Q denotes a generic heavy 
quark, has been described in some detail in Refs. [30,31] (see also Ref. [32]). Following closely 
Refs. [30,31], we have performed two independent calculations of the virtual and real corrections, 
which are in mutual agreement. A detailed account of one of the two calculations of the virtual 
corrections is presented in Ref. [33] . In the following we provide a short summary of our methods 
and mention the tools that have been used. 

The Feynman diagrams and amplitudes that contribute to the virtual corrections have been 
generated with FeynArts 1.0 [34] and FeynArts 3.2 [35]. The amplitudes have been processed 
further with two independent in- house Mathematica routines, which automatically create output 
in Fortran and C++, respectively. The IR (soft and collinear) singularities have been regularized 
in D = 4 — 2e dimensions and have been separated analytically from the finite remainder as 
described in Refs. [31,36]. This separation also allows for a transparent evaluation of rational 
terms that result from D-dependent factors multiplying IR divergences appearing as poles in 
e; in agreement with the general arguments given in Ref. [37] we find that rational terms of 
IR origin cancel completely. The pentagon tensor integrals have been reduced directly to box 
integrals following Ref. [38]. This method does not introduce inverse Gram determinants in the 
reduction process, thereby avoiding numerical instabilities in regions where these determinants 
become small. Box and lower-point integrals have been reduced to scalar integrals using the 
standard Passarino-Veltman technique [39]. Sufficient numerical stability is already achieved in 
this way, but further improvements with the methods of Ref. [40] are in progress. The scalar 
integrals, finally, have been calculated either analytically or using the results of Refs. [41]. The 
IR-finite scalar integrals have furthermore been checked with LoopTools/FF [42]. 

Both evaluations of the real-emission corrections employ (independent implementations of) 
the dipole subtraction formalism [43] for the extraction of IR singularities and for their combina- 
tion with the virtual corrections. Helicity amplitudes for the real emission processes have been 
generated and evaluated with Madgraph [44] and HELAS [45]. The result has been checked by 
an independent calculation using standard trace techniques. 



2.3 Parameter renormalization and resummation improvements 

The renormalization of the strong coupling as(/^) and the factorization of initial-state collinear 
singularities are performed in the MS scheme. As usual, the top quark and the SUSY particles 
are decoupled from the running of as{fJ-)- In the 4FS calculation presented here, also the bottom 
quark is decoupled and the partonic cross section is calculated using a four-flavour Ug. While the 
top- and bottom-quark masses are defined on-shell, the MS scheme is adopted for the renormal- 
ization of the bottom-Higgs Yukawa coupling, which is fixed in terms of the corresponding MS 
renormalization of the bottom mass. In order to sum large logarithmic corrections oc ln(/i/mb) 
we evaluate the Yukawa coupling with the running b-quark mass mb(/i) [46]. 

The SUSY loop corrections induce a modification of the tree-level relation between the bottom 
mass and its Yukawa coupling, which is enhanced at large tan f3 [47-50] . These corrections can 
be summed to all orders by the replacement 

mbtan/3 mbtan/3 (1 — Ab/tan^/3) 

V ^ V (1 + Ab) ^ ^ ^ 
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in the bottom Yukawa coupling [51,52], where 




(2.4) 



with Cp = 4/3 and the auxihary function 



I{a,b,c 



(a2-52)(62-c2)(a2-c2) 



1 



a^b^ In ^ + b^c^ In ^ + c^a^ In . 
0^ / 



(2.5) 



Here, bi 2 are the sbottom mass eigenstates, and mg is the gluino mass. The summation for- 
mahsm can be extended [52] to include corrections proportional to the trilinear coupling A^,. 
However, for the MSSM scenarios under consideration in this work, these corrections turn out 
to be small, and the corresponding summation effects may safely be neglected. 

If the LO cross section is expressed in terms of the bottom Yukawa coupling including the 
summation of the tan /3-enhanced corrections (12. 3p . the corresponding NLO contribution has to be 
subtracted from the one-loop SUSY-QCD calculation to avoid double counting. This subtraction 
is equivalent to an additional finite renormalization of the bottom mass according to 



As we shall demonstrate in the numerical analysis presented in Section El the SUSY-QCD 
radiative corrections are indeed sizeable at large tan (3. After summation of the tan /3-enhanced 
terms, however, the remaining one-loop SUSY-QCD corrections are very small, below the percent 
level. 

3 Phenomenological analysis 

In this section we present NLO SUSY-QCD predictions for the production of heavy charged 
MSSM Higgs bosons at the LHC. We discuss total cross sections and differential distributions 
and compare with the 5FS calculations at NLO for the inclusive tH' cross section. 

3.1 Input parameters 

Let us first specify the values of the input parameters that enter the numerical analysis. 

SM and MSSM masses: The top-quark mass is defined on shell and set to 172.6 GeV [53]. 
For the bottom pole mass we adopt the value mb = 4.6 GeV, corresponding to a MS mass 
mb(?TTb) = 4.26 GeV. The bottom pole mass enters the calculation of the matrix elements 
and the phase space, while the Higgs Yukawa coupling is evaluated using the running 
bottom mass. As for the MSSM parameters, we will focus on the benchmark scenario 
SPS lb [54] which is characterized by a large value of tan/3 = 30 and a correspondingly 
large associated production cross section pp — > tbH^ + X at the LHC. The SPS lb in- 
put parameters are specified in Appendix El The MSSM tree-level relations are used to 
determine the squark masses that enter the SUSY-QCD corrections. The charged-Higgs 
mass is calculated from tan/3 and the mass of the pseudoscalar Higgs, Ma, taking into 
account higher-order corrections up to two loops in the effective potential approach [55,56] 
as included in the program HDECAY [57]. For the Higgs mass determination we use a 
five-flavour Og with as{Mz) = 0.120 [58]. The top quark, the squarks, and the gluino are 
always decoupled from the running of the strong coupling. 




(2.6) 
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Higgs Yukawa coupling: The evaluation of the bottom-Higgs Yukawa coupling, which 
involves the running b-quark mass and the summation of the tan /3-enhanced SUSY-QCD 
corrections through Ab, is also based on a five-flavour Qg with as(Mz) = 0.120. Our default 
choice for the renormalization scale that enters the calculation of the running b-quark mass 
is the average mass of the final-state particles, /u = (m-t -|- mb -|- Mjj-)/3. The scale of Og 
in the summation factor of the Yukawa coupling (cf. Eq. (12. 4p ). on the other hand, is 
determined by the masses of the supersymmetric particles in the loop and is chosen as 
/X = {m^^ + rriy^^ + m^)/3. This scale choice for the effective short-distance contributions 
included in the resummed bottom Yukawa coupling is justified by the recent NNLO results 
for the Ab corrections [59]. 

Hadronic cross section: Our cross-section calculation is defined in the four-flavour 
scheme, i.e. with no b quarks in the initial state. Thus, for a consistent evaluation of 
the hadronic cross sections we adopt the MRST four-flavour parton distribution func- 
tions (pdfs) [60]. The partonic cross section is calculated using the corresponding four- 
flavour as with A*^^) = 0.347 GeV at NLO, except for the Higgs Yukawa coupling which 
is evaluated with a five-flavour Og as explained above. Our default choice for the renor- 
malization and factorization scales that enter the partonic cross section and the pdfs is 
H = (mt -|- mb + Mjj-)/3. Note that the LO cross-section predictions have been obtained 
by using the corresponding LO four-flavour pdf set, a LO with A^^^^ = 0.220 GeV for 
the partonic cross section and a LO running b-quark mass using a LO five-flavour Og with 
as{Mz) = 0.130 [60]. 

3.2 Total cross section and scale dependence 

We first discuss the scale dependence of the total pp — > tbH^ -|- X cross section at the LHC. 
Note that in NLO QCD the cross section for the charge-conjugate process pp — > tbH+ -|- X at 
the LHC is identical to pp tbH~ +X and can be included by multiplying the results presented 
below by a factor of two. The renormalization and factorization scales that enter the hadronic 
cross section and the running b-quark mass are identified and varied around the central value 
/io = (m-t -|- mb -|- -/Vfj£-)/3, but the scale of Og in the summation factor of the Yukawa coupling 
(cf. Eq. (12. 4p ). on the other hand, is kept fixed. Figure [2] (l.h.s.) shows the scale dependence 
of the LO and complete NLO SUSY-QCD cross sections at the SPS lb benchmark point and 
Ma = 200 GeV, corresponding to M^± = 214.28 GeV. As anticipated, the scale dependence 
of the theoretical prediction is significantly reduced at NLO, with a remaining uncertainty of 
approximately ±25% when /x is varied between ^o/3 and 3/xo, compared to approximately ±100% 
at LO. At the central scale, the K-factor, K = crNLo/c"LO is close to one. Note, however, that the 
K-factor strongly depends on the definition of the LO cross section. As described above, our LO 
cross section prediction includes the summation of a certain class of QCD corrections through a 
running Yukawa coupling, and has been evaluated using LO pdfs and ag. The reduction of the 
spurious scale dependence at NLO is particularly striking for the exclusive cross section where 
the b quark is required to be produced with pT,b > 20 GeV, see the r.h.s. of Figure [2l The QCD 
corrections for the exclusive cross section are moderate and negative at the central scale, with a 
corresponding K-factor of K ~ 0.85. 

The total LO and NLO SUSY-QCD cross sections for pp tbH~ +X at the LHC are shown 
on the l.h.s. of Figure[3]as a function of the Higgs-boson mass. Note that tbH~ production at the 
LHC is dominated by gluon-induced processes which provide more than 95% of the cross section. 
The K-factor is displayed in the lower part of the Figure, together with the scale dependence of 
the LO and NLO predictions. We observe that the K-factor is moderate over the whole range 
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2500 




Figure 2: Variation of the LO and NLO cross sections with the renormahzation and factorization 
scales for pp — > tbH~ + X at the LHC, without (l.h.s.) and with (r.h.s.) a cut of pT,b > 20 GeV 
on the b-quark transverse momentum. 




Figure 3: Total LO and NLO cross sections for pp tbH^ + X at the LHC as a function 
of the Higgs-boson mass, without (l.h.s.) and with (r.h.s.) a cut of pT,b > 20 GeV on the b- 
quark transverse momentum. The lower plots show the K- factor, K = crNLo/fLO! ^-nd the scale 
dependence of the LO and NLO cross section predictions for //o/3 < /x < 3^o- 
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Ma [GeV] 


Mh± [GeV] 


mNLO(^) [GeV] 


cj(pp ^ tbH" + X) [fb] 


K — 0-NLo/o"LO 


LO 


NLO 


200 


214.28 


2.80 


583 


562(2) 


0.96 


300 


309.70 


2.76 


248 


251(1) 


1.01 


400 


407.33 


2.72 


114 


119(1) 


1.04 


500 


505.88 


2.68 


56.5 


61.0(2) 


1.09 



Table 1: Total cross sections and K-factors for pp — > tbH~ + X at the LHC. The 
renormahzation and factorization scales are set to = (mt + m-b + Mjj-)/3. The error 
from the Monte Carlo integration on the last digit is given in parenthesis if significant. 
The MRST four-flavour pdfs [60] are adopted. In the third column we show the running 
b-quark mass evaluated at the default renormahzation scale. 



of Higgs-boson masses, with a small increase from K = 0.97 at Mjj± = 200 GeV to K = 1.1 at 
Mjj± = 500 GeV, and that the scale dependence is reduced at NLO also for large Higgs-boson 
masses. Representative values for the total cross section are listed in Table [H To facilitate the 
comparison with other calculations we also show in Table [1] the running b-quark mass, which 
enters the Higgs Yukawa coupling and thus strongly affects the overall normalization of the cross 
section. Requiring the bottom quark to be produced with px,b > 20 GeV reduces the inclusive 
cross section by approximately 60%, see r.h.s. of Figure [31 A systematic comparison of our 
calculation with the results for the exclusive cross section with px,b > 20 GeV and ]%! < 2.5 
presented in Ref. [24] is in progress. 

If we adopt - inconsistently - the five-flavour pdf MRST2004 [58], on which the four-flavour 
set is based, the cross section decreases by approximately 10%: gluon splitting into bottom-quark 
pairs is included in the evolution of the five-flavour pdf and depletes the gluon flux compared to 
the four-flavour pdf. Note that the recent fixed-flavour parton densities of Ref. [61] are based 
on three active flavours in the proton and five active fiavours in the evolution of we can thus 
not use the pdf set of Ref. [61] without modification of our calculation. 

In Table [2] we show the individual contributions to the NLO cross section due to the Standard 
Model QCD corrections and the genuine SUSY-QCD effects, split further into the impact of the 
tan /3-enhanced corrections included in the summation factor Ab and the remainder of the genuine 
SUSY contributions. The cross section labeled ao denotes the LO parton cross section evaluated 
with NLO running b-quark mass, pdfs and Og. The NLO Standard Model QCD corrections, Sqcb, 
increase the prediction by approximately 60%, nearly independent of the value of the Higgs- 
boson mass. This increase is partially compensated by the tan /3-enhanced SUSY corrections, 
'^Sljs'^^'^'^^'^™ ' '^hich amount to approximately —30%. The impact of the remaining one-loop 
SUSY-QCD corrections, ^gug y'^'^'^'^ , is marginal, below the percent level. We also show the result 
of a fixed-order SUSY-QCD calculation, ct^l'o""'^*^'^'^ ' '^liich does not include the tan /3-enhanced 
corrections beyond NLO. We find that the effect of the tan /^-summation beyond NLO, included 
in our best cross-section prediction itnlO; is moderate, at the level of 10%. 

Supersymmetric electroweak 0{a) corrections have been studied in Ref. [62] for the related 
process of neutral MSSM Higgs production in bottom-quark fusion. It has been shown that the 
leading corrections can be taken into account by an appropriate definition of the couplings and the 
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Mh± [GeV] 




fixed— order rrn i 
'^NLO \-^"\ 


TO m 




rtan /3— resum. 
"SUSY 


jrremaindcr 
"SUSY 


214.28 


512 


0.55 


-0.30 


-0.0008 


562(2) 


309.70 


224 


0.61 


-0.30 


-0.0012 


258(1) 


407.33 


106 


0.61 


-0.30 


-0.0009 


125(1) 


505.88 


53.3 


0.62 


-0.30 


-0.0002 


64.1(2) 



Table 2: LO total cross section ctq and NLO corrections 5 relative to do for pp — > tbH~+X 
at the LHC. The error from the Monte Carlo integration on the last digit is given in 
parenthesis if significant. The MRST pdfs [60] are adopted and the renormalization and 
factorization scales have been set to ^ = {nit + "ib + -^h-)/3- "QCD" denotes the 
NLO QCD corrections only, "SUSY/tan/3-resum." the tan /3-enhanced SUSY corrections, 
"SUSY/remainder" the remaining one-loop SUSY corrections and "NLO/fixed-order" the 
complete NLO calculation without summation of the tan /3-enhanced terms. 



running b-mass in an improved Born approximation. The remaining non- universal corrections 
have been found to be small, typically of the order of a few percent. It would be interesting 
to see whether similar conclusions also hold for the process of charged-Higgs production studied 
here. 

3.3 Differential distributions 

Let us now turn to the transverse-momentum and rapidity distributions of the final-state particles 
shown in Figure |H The distributions have been evaluated for the default scale choice fj, = 
(nit + nik + Mh-)/3. The px-distributions of the top quark and the Higgs boson are rather 
similar, with a maximum at ~ 100 GeV. The transverse-momentum distribution of the 
bottom quark is much softer with o"NLo(j'T,b < 25 GeV)/o"NLO ~ 0.7. The heavy particles, 
i.e. the top quark and the Higgs boson, are preferentially produced at central rapidities with 
\y\ ^ 2.5, while the rapidity distribution of the bottom quark is rather fiat in the region \y\ ^ 4. 

The impact of the higher-order corrections on the shape of the transverse-momentum and 
rapidity distributions is shown in Figures [H [H and O respectively. The lower part of the 
Figures shows the K-factor. We find that the shape of the top and Higgs transverse-momentum 
distribution is not strongly affected by the higher-order corrections in the range of j?t relevant for 
the experimental analysis. On the other hand, the bottom quark pT-distribution, which extends 
to pT,h ^ n^h, is softened at NLO, with the K-factor decreasing from K = 1.1 at pT,b ~ 20 GeV 
to K = 0.5 at pT,h ~ 300 GeV. The large impact on the pT,b distribution is due to collinear gluon 
radiation off bottom quarks that is enhanced by a factor Og ln(mb/pT,b)- The enhancement should 
be significantly reduced if the bottom quarks are reconstructed from jets, since the application 
of a jet algorithm treats the bottom-gluon system inclusively in the collinear cone, so that the 
logarithmic enhancement cancels out. The NLO corrections do not significantly change the shape 
of the rapidity distributions at central rapidities \y\ ^ 2, but generically predict more particles 
at large rapidities \y\ ^ 2. 
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Figure 4: NLO transverse-momentum and rapidity distributions of the Higgs boson, the top 
quark, and the bottom quark for pp tbH^ + X at the LHC. 
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Figure 5: LO and NLO transverse-momentum and rapidity distributions of the Higgs boson for 
pp tbH^ + X at the LHC. The lower plot shows the K-factor, K = (Tnlo/clo- 
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100 200 300 400 500 600 -4-3-2-10 1 2 3 4 

PT,t[GeV] yt 

Figure 6: LO and NLO transverse-momentum and rapidity distributions of the top quark for 
pp — > tbH^ -1- X at the LHC. The lower plot shows the K-factor, K = (Tnlo/c'lo- 
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Figure 7: LO and NLO transverse-momentum and rapidity distributions of the bottom quark 
for pp — > tbH~ -|- X at the LHC. The lower plot shows the K-factor, K = cnlo/clo- 
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Figure 8: Total NLO cross section for pp tH~ + X at the LHC as a function of the Higgs- 
boson mass in the 4FS and the 5FS. Shown is the central prediction and the scale dependence 
for ij,q/3 < fi < SfiQ. 

We have also evaluated the differential distributions with the renormalization and factor- 
ization scales set to the average transverse mass, /i = (mT,b + "iT,t + "it,h)/3, where rriT^b = 
Y^'jn.^ + b' ^^'^ ^^^^ shapes of the NLO distributions are not significantly affected 

by such a change. The LO transverse-momentum distributions, however, do provide a better 
description of the NLO shapes when evaluated with fi = (mT,b + ^T,t + "T'T,h)/3. 

3.4 Comparison with the 5FS calculation 

As discussed in Section [H in the 5FS the LO process for the inclusive tH^ cross section is 
gluon-bottom fusion, gb — > tH^. The NLO cross section includes 0{as) corrections to gb 
tH^ and the tree-level processes gg — > tbH^ and qq tbH^, and has been calculated in 
Refs. [21,22]. In Figure [8] we present a comparison of the 4FS and 5FS calculations at NLO 
QCD for the inclusive pp — > tH~ -|- X cross section at the LHC. The 5FS calculation is taken 
from Ref. [21] and is evaluated with the five-flavour pdf MRST2004 [58] and the set of input 
parameters described above. In particular, the renormalization and factorization scales have 
been set to /xq = (mt -|- mb + Mh-)/3, as in the 4FS calculation. The error band indicates 
the theoretical uncertainty when the renormalization and factorization scales are varied between 
/io/3 and S^q. Thus, the error band also includes the scale choice = {rrit + Mjj-)/5 for the 
5FS calculation advocated in Refs. [21,22]. Note that the cross sections shown in Figure [8] do 
not include the NLO SUSY effects, which can be incorporated within good precision by simply 
adjusting the bottom Yukawa coupling according to Eq. (j2.3p . Even taking the scale uncertainty 
into account, the 4FS and 5FS cross sections at NLO are barely consistent; the central predictions 
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in the 5FS are larger than those of the 4FS by approxhnately 40%, rather independent of the 
Higgs-boson mass. 

3.5 Discovery reach 

Accurate theoretical predictions for the charged-Higgs production cross section are crucial to 
exploit the LHC potential for MSSM Higgs-boson searches. To exemplify the importance of 
reducing the theoretical uncertainty through NLO calculations, we consider the discovery reach 
in the search channel pp tbH^ + X followed by the hadronic decay t^v^ with r 

hadrons + as analyzed for the CMS detector in Refs. [25,63]. The number of signal events is 
given by 

iVsignal = //:XC7(pp ^ tbH^ + X) X BR(H^ t^Ut) X BR(r — > hadrons) x exp. efficiency , (3.1) 

where J C denotes the collider luminosity. The experimental efficiency has been determined in 
Ref. [25] as a function of the Higgs-boson mass: 



Mh± [GeV] 


171.6 


180.4 


201.0 


300.9 


400.7 


600.8 


exp. eff. [10-^] 


3.5 


4.9 


5.0 


23 


32 


42 



The QCD background processes lead to 1.7 ± 1 events after cuts, independent of Mjji, so that 
14 or more signal events are needed for a 5 o" discovery [25] . We determine the number of signal 
events from Eq. (j3.ip for the benchmark scenario SPS lb, varying tan/3 and Ma while keeping 
all other supersymmetric parameters fixed. The branching ratio BR(H^ — > t^Vt) varies strongly 
with A/a and has been calculated with SUSY-Hit [64]. The branching ratio of the hadronic r 
decay has been set to BR(r — > hadrons) = 0.65 [65], and we assume an integrated luminosity 
of /£ = 30 fb-^ In Figured we show the 5cr discovery contours for H^ as a function of 
tan/3 and Mjj±, where the number of signal events in Eq. (j3.ip has been evaluated using the 
LO and NLO 4FS calculation presented in this paper. We show results for the central scale 
/io = {Trit + m-b + Mfj-)/3 and results for the renormalization and factorization scales set to 
// = /no/3 and 3;Uo, respectively. Higgs-boson discovery is possible in the areas above the curves 
shown in the figure. Figure O demonstrates that the reduction of the scale uncertainty is crucial 
to exploit the potential of the LHC for charged Higgs-boson discovery. Note that a more detailed 
study of the supersymmetric parameter dependence of the discovery contours is presented in 
Ref. [63] . The importance of a reduced scale dependence through the calculation of higher-order 
corrections for charged-Higgs-boson discovery, however, is generic and largely independent of the 
supersymmetric scenario considered. 

4 Conclusions 

We have presented the next-to-leading order supersymmetric QCD corrections to charged-Higgs- 
boson production at the LHC in the four-flavour scheme through the parton processes qq, gg 
tbH^. While the K-factor is moderate at the central scale /i = (mt -|- rrib + Mh-)/3, the 
QCD corrections considerably reduce the renormalization and factorization scale dependence and 
thus stabilize the theoretical predictions. We find that the shapes of the top-quark and Higgs 
transverse-momentum distributions are not strongly affected by the higher-order corrections. 
On the other hand, the bottom-quark px-distribution is softened at NLO, depending in detail 
on the reconstruction method of the bottom quarks. The NLO corrections do not significantly 
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Figure 9: Discovery reach for MSSM charged Higgs bosons H±, with H± ^ TV, at CMS [25] 
as a function of tan (3 and Mjj± . All other super symmetric parameters have been fixed to the 
SPS lb values. Higgs-boson discovery with J^C = 30 fb~^ is possible in the areas above the 
curves. Shown are results based on the LO and NLO cross sections calculation in the 4FS with 
the central scale ^uq = {rrii + mb + Mj£-)/3 and scales set to ^ = fJ.o/3 and 3fiQ, respectively. 

change the shape of the rapidity distributions at central rapidities, but generically predict more 
particles at large rapidities. We have presented a first comparison of the four-flavour scheme 
NLO inclusive cross sections with a five-flavour scheme calculation based on bottom-gluon fusion. 
Even though the results within the two schemes are consistent within the scale uncertainties, the 
central predictions in the five-flavour scheme are larger than those of the four-flavour scheme by 
approximately 40%. Finally, by referring to a recent CMS study [25] we have demonstrated that 
NLO predictions for the charged-Higgs production cross section are crucial to exploit the LHC 
potential for MSSM Higgs-boson searches. 
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A SPS lb benchmark scenario 



For the SPS lb benchmark [54] scenario discussed in this work we use the following input for 
tan/3, the supersymmetric Higgs mass parameter /x, the electroweak gaugino mass parameters 
Mi^2) the gluino mass nig, the trihnear couphngs AT-^t,h, the scale ;U/j(DR) at which the DR input 
values are defined, the soft SUSY-breaking parameters in the diagonal entries of the squark and 
slepton mass matrices of the first and second generation Mfi (where i = L,R refers to the left- 
and right-handed sfermions, f = q,l to quarks and leptons, and / = n, d, e to up and down 
quarks and electrons, respectively), and the analogous soft SUSY-breaking parameters for the 
third generation M^f: 



tan/3 


= 30.0 




MgL 


= 836.2 


GeV 




= 495.6 


GeV 


MdR 


= 803.9 


GeV 


Ml 


= 162.8 


GeV 




= 807.5 


GeV 


M2 


= 310.9 


GeV 


Mil 


= 334.0 


GeV 


m-g 


= 916.1 


GeV 


MeR 


= 248.3 


GeV 


Ar 


= -195. 


,8 GeV 


M'gE 


= 762.5 


GeV 


At 


= -729. 


,3 GeV 


ME 


= 780.3 


GeV 


A^, 


= -987. 


,4 GeV 


Ml% 


= 670.7 


GeV 


M/j(DR) 


= 706.9 


GeV 


Mf^ 


= 323.8 


GeV 




= 218.6 


GeV. 









The mass of the CP-odd Higgs boson Ma is varied and taken as input to calculate the charged- 
Higgs boson mass Mjji, taking into account higher-order corrections up to two loops in the 
effective potential approach [55,56] as included in the program HDECAY [57]. 
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